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We propose a version of a gravitational theory with the torsion field, induced by the chameleon 
field. Following Hojman et al. Phys. Rev. D 17, 3141 (1976) the results, obtained in Phys. Rev. D 
90, 045040 (2014), are generalised by extending the Einstein gravity to the Einstein-Cartan gravity 
with the torsion field as a gradient of the chameleon field through a modification of local gauge 
invariance of minimal coupling in the Weinberg-Salam electroweak model. The contributions of the 
chameleon (torsion) field to the observables of electromagnetic and weak processes are calculated. 
Since in our approach the chameleon-photon coupling constant is equal to the chameleon- 
matter coupling constant /?, i.e. /I 7 = Pj the experimental constraints on /?, obtained in terrestrial 
laboratories by T. Jenke et al. (Phys. Rev. Lett. 112, 115105 (2014)) and by H. Lemmel et al. 
(Phys. Lett. B 743, 310 (2015)), can be used for the analysis of astrophysical sources of chameleons, 
proposed by C. Burrage et al. (Phys. Rev. D 79, 044028 (2009)), A.-Ch. Davis et al. (Phys. Rev. 
D 80, 064016 (2009) and in references therein, where chameleons induce photons because of direct 
chameleon-photon transitions in the magnetic fields. 

PACS numbers: 03.65.Pm, 04.62.+v, 13.15.+g, 23.40.Bw 


I. INTRODUCTION 


The chameleon field, the properties of which are analogous to a quintessence 0,0 , i.e. a canonical scalar field 
invented to explain the late-time acceleration of the Universe expansion has been proposed in In order 

to avoid the problem of violation of the equivalence principle 0 a chameleon mass depends on a mass density p of a 
local environment 0^ ■ The self-interaction of the chameleon field and its interaction to a local environment with a 
mass density p are described by the effective potential Ves{4>) 


V,s{cj^) = V{cj^)+pe^*/^^\ 


( 1 ) 


where ()) is a chameleon field, d is a chameleon-matter field coupling constant and Mp\ = I/v^SttGjv = 2.435 x 10^^ eV 
is the reduced Planck mass [1^. The potential defines self-interaction of a chameleon field. 

As has been pointed out in Ref. ultracold neutrons (UCNs), bouncing in the gravitational field of the Earth 

above a mirror and between two mirrors, can be a good laboratory for testing of a chameleon-matter field interaction. 
Using the solutions of equations of motion for a chameleon field, confined between two mirrors, there has been found 
the upper limit for the coupling constant f3 < 5.8 x 10® [l^, which was estimated from the contribution of a chameleon 
field to the transition frequencies of the quantum gravitational states of UCNs, bounds in the gravitational field 
of the Earth. For the analysis of the chameleon-matter field interactions in R,efs. (ltl - ll^ the potential VJ6) of a 
chameleon-field self-interaction has been taken in the form of the Ratra-Peebles potential (see also 0 Q) 


V{(^) = A^ + 


A4+r, 


( 2 ) 


where A = = 2.24(2) x 10"® eV [T^ with^ = 0.685p;!^i^ and Hq = 1.437(26) x 10"®® eV are the 

relative dark energy density and the Hubble constant [1^, respectively, and n is the Ratra-Peebles index. The 
runaway form A^+”/(/)” for (p —>■ 00 is required by the quintessence models 0i- Such a potential of a self-interaction 
of the chameleon field allows to realise the regime of the strong chameleon-matter coupling constant /3 ^ 10® (Tfll - [l3| . 

Recently [l^ some new chameleon-matter field interactions have been derived from the non-relativistic approx¬ 
imation of the Dirac equation for slow fermions, moving in spacetimes with a static metric, caused by the weak 
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gravitational field of the Earth and a chameleon field. The derivation of the non-relativistic Hamilton operator of the 
Dirac equation has been carried out by using the standard Foldy-Wouthuysen (SEW) transformation. There has been 
also shown that the chameleon field can serve as a source of a torsion held and torsion-matter interactions [I3-[2l. 
A relativistic covariant torsion-neutron interaction has been found in the following form [l^ 


^t{x) = ^grTf,{x)ij{x)a'"''K'fp{x), ( 3 ) 

^ _ y 

where 7^ is the torsion held, i’ix) is the neutron held operator, A{x)di, B{x) = A{x)duB{x) — {d^A{x))B{x) and 
is one of the Dirac matrices [1^. In the non-relativistic limit we get 

5Cr{x) = igrT■ '!/j(a:)(S x V )tj}{x) + ... = igrT ■ (p\x){a x V )(p{x) + ... , (4) 

where <p{x) is the operator of the large component of the Dirac bispinor held operator '0(x) and E = y°yy^ is the 
diagonal Dirac matrix with elements S = diag((j, tr) and a are the 2x2 Pauli matrices [2^. As has been found in 
[l^ the product g-rT is equal to 

where 7 = 1 for the Schwarzschild metric of a weak gravitational held [l^( see also [26l - [^ l. 

Following [ 1 ^ we introduce the torsion held tensor as follows 

- <5“^/,., ( 6 ) 

where = d^f and / = e'^“. Such an expression one obtains from the requirement of local gauge invariance of the 
electromagnetic held strength [l^ (see also section Hill) . According to [l^, the scalar held (/)h can be identihed with 
a chameleon held ^ as (/)h = /3</'/AIpi. As a result, we get 


rOL _ 

— 


M, 


PI 




/3 


M, 


PI 






(7) 


where = g'^^gxu and g,j\ and g^°‘ are the metric and inverse metric tensor, respectively. The torsion tensor held 
is anti-symmetric 

For the subsequent analysis we need a dehnition of the covariant derivative V^-i, of a vector held in the curve 
spacetime. It is given by [29l - [3ll| 


Vf,,, = ( 8 ) 

where r“^;, is the affine connection, determined by [3 

= {%.} - i - V. - = {%.} + ^ + g^.M. (9) 

where fiv} are the Christoffel symbols [2§-[3l| 

{ = -^g {g\n,v + g\u,ti ~ gtj,v,\) (lO) 

and is the torsion tensor held [Bill 




T a _ -pQ: _ po: _ 


PI 




( 11 ) 


We introduce the contribution of the torsion held in agreement with Hojman et al. [l^ (see Eq.(38) of Ref. 0 ). 

Having determined the affine connection we may introduce the Riemann-Christoflel tensor TZA or curvature 
tensor as [29l - [3lj| 


‘T)Ct -pQ; -pQ! I p'Q; p't^ pia pii^ /i o^ 

^\,y \ ^ fiiy -L fiXi 

which is necessary for the dehnition of the Lagrangian of th e g ravitational held in terms of the scalar curvature 77 
[ 2 ^ related to the Riemann-Christoffel tensor TZA^vX by [29l - [31j 


n = g>^^n‘^^^x = g^^n^x- 


( 13 ) 
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Here Ti.^\ is the Ricci tensor [2914^ . Following [T^ and skipping intermediate calculations one may show that the 
scalar curvature TZ is equal to 

n = (14) 

where the curvature R is determined by the Riemann-Christoffel tensor Eq. (HU with the replacement ^ {“/ii/}- 

The paper is organised as follows. In section |lT] we consider the chameleon field in the gravitational field with 
torsion (a version of the Einstein-Cartan gravity), caused by the chameleon field. We derive the effective Lagrangian 
and the equations of motion of the chameleon field coupled to the gravitational field (a version of the Einstein gravity 
with a scalar self-interacting field). In section Hill we analyse the interaction of the chameleon (torsion) field with the 
electromagnetic field, coupled also to the gravitational field. Following Hojman et al. [I^ and modifying local gauge 
invariance of the electromagnetic strength tensor field we derive the torsion field tensor T°‘ in terms of the chameleon 
field (see Eq. ([T]). In section IIVI we analyse the torsion (chameleon) - photon interactions in terms of the two-photon 
decay ^ 7 + 7 of the chameleon and the photon-chameleon scattering 7 + </> —>• </> + 7 . We show that the amplitudes 

of the two-photon decay and the photon-chameleon scattering are gauge invariant. In order words we show that the 
replacement of the photon polarisation vectors by their 4-momenta leads to the vanishing of the amplitudes of the 
two-photon decay and the photon-chameleon scattering. In section|V]'we investigate the Weinberg-Salam electroweak 
model [ 1 ^ without fermions. We derive the effective Lagrangian of the electroweak bosons, the electromagnetic field 
and the Higgs boson coupled to the gravitational and chameleon field. Such a derivation we carry out by means 
of a modification of local gauge invariance. In section IVII we include fermions into the Weinberg-Salam model and 
derive the effective interactions of the electroweak bosons, the Higgs held and fermions with the gravitational and 
chameleon Held. In section Eni we calculate the contributions of the chameleon to the charge radii of the neutron 
and proton. We calculate the contributions of the chameleon to the correlation coefficients of the neutron /3“-decay 
n—>p + e“ + i>e + <(> with a polarised neutron and unpolarised proton and electron. In addition we calculate the cross 
section for the neutron /3“-decay ()) + n —>■ p + e“ + 1/e, induced by the chameleon Held. In section rVlIIl we discuss 
the obtained results and perspectives of the experimental analysis of the approach, developed in this paper, and of 
observation of the neutron /3“-decay, induced by the chameleon. 


II. TORSION GRAVITY AND EFFECTIVE LAGRANGIAN OF CHAMELEON FIELD 

The action of the gravitational Held with torsion, the chameleon Held and matter Helds we deHne by 0 

5'g.ch = J + j d!^x^/^Cnr[g^,y\, (15) 

where the Lagrangian C[R,(j)\ is given by 

C[n, 0] = i Mil i (1 - 3 / 32 ) _ V{4>). (16) 

Here = d(j)/dx^ and (j)’^ = d(i)/dx^ and V{4i) is the potential of the self-interaction of the chameleon Held Eq.©. 
The matter Helds are described by the Lagrangian Cm[g^^v\■ The interaction of the matter Held with the chameleon 
Held runs through the metric tensor in the Jordan-frame i a m, m, which is conformally related to the 
Einstein-frame metric tensor by g^^ = P g^v (or g^^'^ = f~'^ g^'') and with / = p. [^. 

The factor can be interpreted also as a conformal coupling to matter Helds [1,[3. Using Ea. dHl) we transcribe 

the action Eq. (HU into the form 

S'g.ch = J Pxyf^ (^^MpiR +C[(l)]^ + j Pxy/^CmlgtJ.,^], (17) 

where the contribution of the torsion Held to the scalar curvature is absorbed by the kinetic term of the chameleon 
Held. The Lagrangian £[</>] is equal to 

m = (18) 

The total Lagrangian in the action Eq. m is usually referred as the Lagrangian in the Einstein frame, where 
as well as g^" is the Einstein-frame metric such as = g^'^(j )and □</> = {i/y/^){y/^4>’P-,ti = 

[11 (see also [i,[3). 
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Varying the action Eq. (13 with respect to (j)^^ and (j) we arrive at the equation of motion of the chameleon field 

d S{^/^C[(j)]) ^ S{y/^C[(I)]) 5{yf^Cm) 

dxf^ 64> 6cj) 

Using the Lagrangian Eq. m we transform Eq. (m into the form 


d 




2 S(y/^£m) 


where V^{(j)) and are derivatives with respect to (/). Since by definition [1,[3 the derivative 

2 - 

is a matter stress-energy tensor in the Jordan frame, Ea. (Eni) takes the form 


( 20 ) 


( 21 ) 


( 22 ) 


where T“ = g‘^^Ta\. For a pressureless matter T“ = ^where p is a matter density in the Jordan frame, related to 
a matter density in the Einstein frame p hy p = f~^ p [6| we get 

□0 = -U;(0)-p/;, (23) 

where we have set = P 0- Then, V^{(j)) — p fIp coincides with the derivative of the effective potential of the 

chameleon-matter interaction Ueff(</>) with respect to (p, given by Eq.(IT|) for / = 


III. TORSION GRAVITY WITH CHAMELEON AND ELECTROMAGNETIC FIELDS 


In this section we analyse the interactions of the torsion (chameleon) field with the electromagnetic field. The 
action of the gravitational field, the chameleon field, the matter fields and the electromagnetic field is equal to 


5, 


g,ch,em 


+ 


d X \l g £m ’ 


(24) 


where = p Since = p we get that = P^g°‘^g^''- The term P^CPg^v] describes 

an environment where the chameleon field couples to the electromagnetic field. 

Following then Hojman et al. [l^ we define the electromagnetic strength tensor field in the gravitational and 
torsion field 




fll' 


— A —A —A —A —A 


- A^T 


fiiy 


1 


(25) 


where — A^^i, and A^ is the electromagnetic 4-potential. According to Hojman et al. [I^, under a gauge 

transformation the electromagnetic potential transforms as follows 


Afj, —>■ A^ — Af^ + c^“($)A_a, ( 26 ) 

where c^“($) is a functional of the scalar field 4>, which we identify with the chameleon field 4> = (3(j)/Mp\, i.e. 
c^“($) cpp) = and A is an arbitrary gauge function. The gauge invariance of the electromagnetic 

field strength imposes the constraint [T^ 


cPP),^ - cpp),, - cPp)r\. = 0. (27) 

This gives the torsion tensor field given by Eq.® and Eq.®. Substituting Eg. (1^51) into Ea. ip^ we arrive at the 
expression 

5g.ch.em = JdPP^ (i i? + i - Vp) - i 


(28) 
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Using the definition of the torsion tensor field Eq.® we arrive at the action 

5g,ch,em = Jd^x^ (i M|ii? + i - U((^) - i - AM 

- \^2^A^^’''-A''4>'''){A^,(i)^^,-A^(j)^^,)^+jd'^xM^F.m[gt,u] ( 29 ) 

Thus, because of the torsion-electromagnetic field interaction the chameleon becomes unstable under the two-photon 
decay ()>—>■ 7 + 7 and may scatter by photons 7 + ^ —?►(/) + 7 with the chameleon-matter coupling constant /3/Mpi. 
These reactions are described by the effective Lagrangians 

M -ir ( 30 ) 

Mpi 


and 


1 /3^ 
.1 


- AM = 




( 31 ) 


For the application of the action Eq. (ESI) with chameleon-photon interactions to the calculation of the specific reactions 
of the chameleon-photon and chameleon-photon-matter interactions we have to fix the gauge of the electromagnetic 
field. We may do this in a standard way 


5g.ch.em = Jd^xM Mil i? + 1 cj^M - y{4>) - \ FM" + fM 

+ J d'^x C 

m 

where f {A^^, A^, (f)) is a gauge fixing functional and the divergence A>^.^ is defined by Him 

\T^d.AV 


We get 


The afhne connection f we have to calculate for the Jordan-frame metric 

f“. = {%.} + <5% /-' /.. = {%.} + ^ s% 

such as = T“^j/ (see Eq.d?])). As a result, the divergence is equal to (Min 

^/^ dxf^ Mpi 

Since a gauge condition should not depend on the chameleon field, we propose to fix a gauge as follows 




( 32 ) 


( 33 ) 


( 34 ) 


( 35 ) 


( 36 ) 


where ^ is a gauge parameter. Now we are able to investigate some specific processes of chameleon-photon interactions. 


IV. CHAMELEON-PHOTON INTERACTIONS 

The specific processes of the chameleon-photon interaction, which we analyse in this section, are i) the two-photon 
decay (j) —>■ 7+7 and ii) the photon-chameleon scattering j+(j) —>■ ^(>+ 7 . The calculation of these reactions we carry out 
in the Minkowski spacetime. For this aim in the interactions Eq. (111(11) and Eq. (HTTI) we make a replacement —>■ , 

where g^''' is the metric tensor in the Minkowski space time with only diagonal components (+1, —1, —1, —1), and 
g -)■ det{77i,j.} = -1. 
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FIG. 1: Feynman diagram for the 7 + 7 decay. 
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FIG. 2: Feynman diagram for the amplitude of the photon-chameleon (Compton) scattering 


A. Two—photon 


7 + 7 decay of the chameleon 


For the calculation of the two-photon decay rate of the chameleon we use the Lagrangian Eg. (1301) . The Feynman 
diagram of the amplitude of the two-photon decay of the chameleon is shown in Fig.[TJ The analytical expression of 
the amplitude of the (/>—?► 7 + 7 decay is equal to 


M(<^ ^ 77 ) = - 2 ^ ((£* . e\){k^ • fc 2 ) - • k^){el • fci)), 

Mpi 


(37) 


where e*{kj) and kj for j = 1,2 are the polarisation vectors and the momenta of the decay photons, obeying the 
constraints E*{kj) ■ kj =0. Skipping standard calculations we obtain the following expression for the two-photon 
decay rate of the chameleon 


ml 


where is the chameleon mass, defined by [TTj | 


n+2 

2n+2 


(38) 


(39) 


as a function of the chameleon-matter coupling constant /3, the environment density p and the Ratra-Peebles index 
n. 

B. Photon—chameleon 7 + (()—><() + 7 scattering 

The Feynman diagrams of the amplitude of the photon-chameleon scattering are shown in Fig. [2] The contributions 
of the diagrams in Fig. [5] are given by 


M(“)(7()) ())7) = ^ (k>^{ei-pi)-e'({ki-pi)^ Df,c.{q)(k^{el ■ P 2 ) - el°^{k2 ■ P2)j 

+ ^ - q>^q^{q) - q'^q<^D^P{q) + q’'q^(q)) e*^^P2p 

+ (fci(ei - Pi) -ei(fci - Pi)) {q“D/(q) - q^Df,°‘{q)^ 


9=Pl+fcl=P2+fe2 

9=Pl+fcl=P2+fe2 


/3^ 




(^ 2(^2 ■P 2 ) - ETik 2 -^ 2 )) 


g=pi+fei=P2+fe2 


5=pi+fcl=P2 + fc2 


(40) 
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(^2(^2 -Pi) - -Pi)) Df,aiq)(k“{ei ■P2) - e^{ki -^2)) 


M 2 

0 


q=pi-k2=P2-ki 

+ ^ e*^^p,.(^q>^q‘^D‘'^iq) - q>^q^D'^^{q) - q^q^{q) + q^q^{q)) e,^P2p 
^ (^ 2(^2 - Pi) - £^(*2 - Pi)) {q°‘Df{q) - q^D 0 ^{q)^ eiaP 2 P 


q=Pl-k2=P2-ki 


M 2 , 


M 2 , 


4 pPi<^{q^D'"a{q) - q''D^"c^ (^f(ei ^2) - e?(A:i •P2)) 


<jr=pi-fc 2 =p 2 -fcl 

q=pi-k2=P2-ki 


(41) 


m(°)(7(() ()!)7) = - 2 (e; • ei){pi ■ P2) + (£2 •Pi)(£i •P2) + (£2 •P2)(£i -Pi)), 


(42) 


= 2 n(n + l)(n + 2 ) ^ ^ 

^P1 CTn m^-q^- ^0 


q=k2-kx=pi-p2 


(43) 


where £1 and £2 are the photon polarisation vectors in the initial and final states of the photon-chameleon scattering. 
They depend on the photon momenta £i(fci) and £ 2 (^ 2 ) and obey the constraints £i(A:i) • ki = £ 2 (^ 2 ) • ^2 = 0. The 
chameleon field mass is defined by Eg . (1391) . The vertex of (f0 interaction is defined by the effective Lagrangian 




i(n+ l)(n + 2 ) A"+4 


^n+3 

min 


Here 4>-mm is the minimum of the chameleon field, given by Eini 

, . / ?^Mpi\ ;7+T 


(44) 


(45) 


where p is the density of the medium in which the chameleon field propagates. The photon propagator Dap{q) is 
equal to 


Dap{q) — ^2 


to 


(«-(i-a 2 ^). 


( 46 ) 


One may show that the amplitudes (j) —>■(/> 7 ) and cj) —>■(/> 7 ) do not depend on the longitudinal part of 

the photon propagator. As a result the amplitudes —>■(/) 7 ) and (j) —^ 7 ) can be transcribed into the 

form 

^ ((£^P 2 )(£i ■pi)(A:i ■ ^ 2 ) - (£^P 2 )(£i • fc 2 )(fci-Pi) 


/32 

1 

Ml, 

q^ + *0 

-{e *2 

• fcl)(£l 

/32 

1 

M|, 

q^ + iO 

-(£2 

■ 9 )(£i ■ 

/32 

1 

M|i 

q^ + to 

-(£; 

• fcl)(£l 

0 ^ 

1 

M|, 

q^ + iO 

-{^2 

■Pi)(£i 


(J=Pl+fel=P2+fc2 


(J=pi+fcl=P2+fc2 


((£2 • 9 )(£i • Pi)(fei ■ P2) - (£2 • '?)(£! ■ P2)(fci • Pi) 


g=pi+fci=P2+fe2 


((£2 ■P2)(£i ■ q){k2-pi) - (£2 ■P2)(£i • k2){pi ■ q) 


-(£2 ■Pi)(£i ■ q){k2 •P2) + (£2 • £i)(fc2 •P2)(P1 • 9)) 


q=pi+fcl=P2+fe2 


(47) 
































and 


(j)j) = ^ ^ (^(S2 ■ Pl)(si ■ P2)(ki ■ fcz) - (£2 ■ fcl)(£l ■ P 2 )(k 2 ■ Pi) 


+ 


-(£2 • Pi)(£i • k2)iki ■ P2) + (£2 • ei){ki ■ P 2 ){k 2 ■ Pi)) 

((£^2 ■ 9 )(ei ■q){Pi ■P 2 ) - (£2 ■ q){£i ■ Pi){P 2 ■ q) 


/32 

1 

M|i 

g2 + iO 

-(£2 

■Pi){.£i 

/32 

1 

M|, 

g2 + iO 


■P 2 ){£i 

/32 

1 

M|i 

g2 + iO 

-(^2 

■Pi){£i 

/32 

1 

Mil 

g2 + iO 

-(^2 

■q)[£i ■ 


<jr=pi-fc2=P2-fel 


q=pi-k2=P2 — kl 


^ ((i£2 • Pl)(.£l ■ q){k2 ■ P2) - (^2 • P2){,£l ■ q){k2 ■ Pi) 


q=Pl-k2 =P 2 -ki 


-(£2 • q){£i ■ Pi)iki ■ P2) + (£2 • £i){ki ■ P 2 )ipi ■ <?)) 


q=Pl—k2 =P 2 -ki 

The total amplitude of the photon-chameleon scattering is defined by the sum of the amplitudes Eg. lHO)) 


(48) 


M{'j (j) ^ (j)'y) = (j) ^ (j)'y). (49) 

j=a,b,c,d 

Now let us check gauge invariance of the amplitude of the photon-chameleon scattering Ea. lH^ . As we have found 
already the amplitudes M^°'\'y (j) —>•()) 7 ) and M^^\'y (j) —>■ <^ 17 ) do not depend on the longitudinal part of the photon 
propagator, i.e. on the gauge parameter Then, according to general theory of photon-particle ( 7 / 1 ) interactions 
[ 23 , the amplitude of photon-particle scattering should vanish, when the polarisation vector of the photon either in 
the initial or in the final state is replaced by the photon momentum. This means that replacing either £1 —>■ ki or 
£2 —>■ ^2 one has to get zero for the total amplitude Ea. (H^ . Since one may see that the amplitude M'^'^\'y (j) —>■ ^iy) is 
self-gauge invariant, one has to check the vanishing of the sum of the amplitudes, defined by the first three Feynman 
diagrams in Fig.[2l i.e. 


Replacing £1 —>• fci we obtain 
M^°'\'y 4> —>■(() 7 ) 
M^^\'y(f> —>■ ^ 7 ) 
M^‘^\'y (j) —>■ ^ 7 ) 
Because of the relation 


£1 — 


£l->fel 


£l->fel 


M (7 ())—>■()) 7) = M^^\'y (j) ^ (j)'y). 

j—a,b,c 


( - P 2 ){ki ■ pi) + {e *2 ■ fci) 


M 2 , 


Mil 


M|i 


q=pi+ki=p2+k2 

q=Pl-k2 =P 2 -ki 


( - Pi){ki ■ P2) + {£*2 ■ ki) 

( - 2(£2 • ki)(pi ■ P2) + (£2 • P2){ki ■ Pi) + (£2 • Pi){ki ■ P2)) • 


5=pi+fcl=P2 + fc 2 

the sum of the amplitudes Ea. dCTl vanishes, i.e. 


+ q^ 


q=pi—k2=P2-ki 


= 2 (pi ■P2) 


M^-^^ (7 ())—>■ ^7) 


j—a^b,c 


ei-)-fei 


= 0 . 


(50) 


(51) 

(52) 

(53) 


The same result one may obtain replacing £2 —^ k 2 - Thus, the obtained results confirm gauge invariance of the 
amplitude of the photon-chameleon scattering, the complete set of Feynman diagrams of which is shown in Fig. [2] 
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Of course, because of the smallness of the constant < 10“®° barn/eV^, estimated for /3 < 5.8 x 10® [T^, the 

cross section for the photon-chameleon scattering is extremely small and hardly plays any important cosmolo gica l 
role at low energies, for example, for a formation of the cosmological microwave background and so on |55l . I51 tI |. 
Nevertheless, the observed gauge invariance of the amplitude of the photon-chameleon scattering is important for 
the subsequent extension of the minimal coupling inclusion of a torsion field to the Weinberg-Salam electroweak 
model [ 2 ^ in the Einstein-Cartan gravity. One of the interesting consequences of the observed gauge invariance of 
the chameleon-photon interaction might be unrenormalisability of the coupling constant /3/Mpi by the contributions 
of all possible interactions. This might mean that the upper bound on the chameleon-matter coupling constant 
/3 < 5.8 X 10®, measured in the qBounce experiments with ultracold neutrons [l^ . should not be change by taking 
into account the contributions of some other possible interactions. 

In this connection the results, obtained in this section, can be of interest with respect to the analysis of the contribu¬ 
tions of the photon-chameleon direct transitions in the magnetic field to the cosmological microwave background . 
The effective chameleon-photon coupling constant ^eff, introduced by Davis, Schelpe and Shaw [s^, in our approach 
is equal to ^eff = PjMpi. Using the experimental upper bound /3 < 5.8 x 10® we obtain geff < 2.4 x 10“^° GeV“^. This 
constraint is in qualitative agreement with the results, obtained by Davis, Schelpe and Shaw The experimental 
constraints on the chameleon-matter coupling j3 < 1.9 x 10^ (n = U, ^5 < 5.8 x 10^ (n = 2), f3 < 2.0 x 10® (n = 3) 
and /3 < 4.8 x 10® (n = 4), measured recently by H. Lemmel et al. |57l | usin g th e neutron interferometer, place more 
strict constraints of the astrophysical sources of chameleons, investigated in |5l|-(^. 


V. TORSION GRAVITY AND WEINBERG-SALAM ELECTROWEAK MODEL WITHOUT 

FERMIONS 


In this section we investigate the Weinberg-Salam electroweak model without fermions in the minimal coupling 
approach to the torsion field (see (T^); caused by the chameleon field [l^ . 

According to Hojman et al. [l^ . in the Einstein-Cartan gravity with a torsion field, induced by a scalar field, the 
covariant derivative of a charged (pseudo)scalar particle with electric charge q should be equal to 


Dtj.=df,-iqf ^ Afj 


(54) 


with / = Using such a definition of the covariant derivative one may calculate the electromagnetic held 


strength tensor as follows 




= -^/[D^,D,] = /(a^(/-iA,) 


— A 


(55) 


where the torsion tensor Held 7)j“ is given by Eg. lfTTl) . In this section we discuss the Weinberg-Salam electroweak 
model psj l in the Einstein-Cartan gravity with a torsion Held, caused by the chameleon Held. Below we consider the 
Weinberg-Salam electroweak model without fermions. 

The Lagrangian of the Weinberg-Salam electroweak model of the electroweak bosons and the Higgs Held with gauge 
SU{2) X 1/(1) symmetry, determined in the Minkowski space-time, takes the form [ 2 ^ 


r ——It . Tmi' _ i h 

-^ew — ^ ^ 

1 


i g' Yu, i 


X - i -g" - i g - U($1'4>), 


(56) 


where g' and g are the electroweak coupling constants and and are vector Helds and $ is the Higgs boson Held. 
Then, and Iw = h Yu are the weak hypercharge and the weak isopin, respectively: fu, = (r^, are the weak 


isospin 2x2 Pauli matrices + i and tr(r“r.^) = 2 [^. The weak hypercharge Y^u and the third 

component of the weak isospin 1^3 are related by Q = + Yu,/2, where Q is the electric charge of the Held in units 

of the proton electric charge e [25| . In the Weinberg-Salam electroweak model the Higgs boson Held $ possesses the 
weak isospin /„, = 1/2 and the weak hypercharge Y^, = 1. The field strength tensors A^i, and B^i, are equal to 


Aful/ - d^Ay OyAf^ g Afu X Ay^ 


Bfuy = df,By - dyBfj,. 


(57) 


The Higgs boson field and its vacuum expectation value are given in the standard form 

$+ \ 1 


$ = 






(58) 
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where = (t; + (p)l\/2 and (/? is a physical Higgs boson field. The potential energy density has also the 

standard form [2^ 

+ (59) 

with > 0, K > 0 and = y? jn. The vacuum expectation value is related to the Fermi coupling constant Gp 
Gpv^ = 1, where Gp = 1.16637(1) x 10”^^MeV“^ [l^. The covariant derivative of the Higgs field is given by 

D^ = d^,-i]^g'B^<^-i]^gT- (60) 
Using the covariant derivative Eg. (1601) we may calculate the commutator [7?^, Di^] and obtain the following expression 
[D^,D^] = -i^g'{d^B^ - d„B^) -ig^r- (5^X - d^Af,+gAf, x A^) = -i^g' B^,y - i^gr- A^^, (61) 


where and are the field strength tensors Eg. (1571) . Under gauge transformations 

A^^^A^ = QA^n-^ + -d^nn-\ 

'^9 

B^^^B^ = B^+d^A, (62) 


where U and A are the gauge matrix and gauge function, respectively, and 
Afiv = Afj^i, and transform as follows 


Bf_ii, 




nAfj,M \ 

Bfiv 


■ Af^, the field strength tensors 


(63) 


In the Einstein-Cartan gravity with a torsion field in the minimal coupling approach the covariant derivative Eo. (l60l) 
should be taken in the following form 


D^ = df,-i-g' f ^ B^-i-gf- 


•• A 




(64) 


Eor the definition of field strength tensors A^j^ and H^i/, extended by the contribution of a torsion field, we propose 
to calculate the commutator [Dfj_,D^]. The result of the calculation is 

[Df,, D^] = -i i g' /-I B^^-i^gf~^T- (65) 

where the field strength tensors and are equal to 

Afj^i/ = Ai/ p^ 'A g f A^ X Ai/ AqAT 

( 66 ) 

where the torsion tensor field fiv is given in Eq. (CD- Thus, the Lagrangian of electroweak interactions in the 
Einstein-Cartan gravity with a torsion field in the minimal coupling constant approach and the chameleon field, 
coupled through the Jordan metric g^n = g^iv, takes the form (34j 

^ + f(d^‘^-i\g’r^B^^-^]^gr^T■A^‘^)^ 

X - i\g' - i]^g T ■ A^^^) - (67) 

where the factor comes from \/—g = y/—g- The physical vector boson states of the Weinberg-Salam electroweak 

model are [l^ 


IU± = -^{AItiAD, 

= sin 9w Bfj^ — cos 6\y A^, 
A^ = cos 9w B^ + sin 9w A^, 


(68) 
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where and are the electroweak W^-boson and Z-boson fields and is the electromagnetic field, respectively, 
and Ow is the Weinberg angle defined by tandw = g'/g. The electromagnetic coupling constant e as a function of the 
coupling constants g and g' is given by e = gg' j \Jg^ + g"^ = g sin^w = g' cos 6w 113 • In terms of the electroweak 
boson fields and the electromagnetic field the Higgs boson field (p and the chameleon field (j), coupled to 
the gravitational field with torsion, the Lagrangian of the electroweak interactions takes the following form 


4 ^ = j i W ++ I Ml W+W->^ + ImI 

V~9 4 4 


sinf 


2 ■ ’ f-iu' * '2 M '2 

Ow - cos6»w {w-'^Z'' - Z^^W-'^ 

sin9w {w^^A'' - A^^W+'''^ - cos9w (w+^Z'' - Z^^W+'^') 

f-^tg(^sm9wA^,- cos 9w {w->^W+'' - W 
f-^g^ (W-Wl - W-Wl^ {W+>^W-'^ - W+^W -4 
r^g'^ [sin^w (W^A, - - cos^h^ (wlZ, - 

sm9w (w-i^A^ - Af^W-''^ - cos9w 


a\ 9M^P WlW-^ + i , i 


^ MzpZf,Zf^+^ 


9 


8 cos^ 9w 




+ \f - MIp^ - f KVP^ - K 


(69) 


where Ml = g^v^ j A:, M| = Ml/cos^9w and M'^ = 2 ku^ are the squared masses of the W^-boson, Z-boson and 
Higgs boson field, respectively, A^iu, and are the strength field tensors of the electromagnetic, Z-boson and 
IT^-boson fields. They are equal to 


A — A —A —A 

= Z,y_^ — Z^^n — Z^T^ 

y\;± = 


Now we are able to extend the obtained results to fermions. 


(70) 


VI. TORSION GRAVITY WITH CHAMELEON FIELD AND WEINBERG-SALAM ELECTROWEAK 

MODEL WITH FERMIONS 

A. Dirac fermions with mass m in the Einstein—Cartan gravity, coupled to the chameleon field through the 

Jordan metric g^v 

The Dirac equation in an arbitrary (world) coordinate system is specified by the metric tensor g^lx). It defines 
an infinitesimal squared interval between two events 

ds^ = g^i,{x)dx^dx''. (71) 

The relativistic invariant form of the Dirac equation in an arbitrary coordinate system is 

(i7^(a:)V^ - m)il}{x) = 0, (72) 

where 7 ^(x) are a set of Dirac matrices satisfying the anticommutation relation 

7'^(x)7"(x) + 7"(x)7'^(x) = 2g^"(x) (73) 

and V is a covariant derivative without gauge fields. For an exact definition of the Dirac matrices 7 ^(a:) and the 
covariant derivative we follow [ 2 ^ and use a set of tetrad (or vierbein) fields e“(x) at each spacetime point x 
defined by 


da;“ = e“(a:)dx^. 


(74) 
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The tetrad fields relate in an arbitrary (world) coordinate system a spacetime point x, which is characterised by the 
index /r = 0,1, 2,3, to a locally Minkowskian coordinate system erected at a spacetime point x, which is characterised 
by the index d = 0,1, 2, 3. The tetrad fields e“(a:) are related to the metric tensor g^j.v{x) as follows: 


dx°‘d.x^ = Va0 [e'^ix)dxf"][e^{x)dx''] = [ri^ge‘^{x)e^{x)]dx^dx’' = g^^{x)dxf^dx''. 


I&0 


&0 


This gives 


= g&p e“(a;)e^(a;). 


(75) 


(76) 


Thus, the tetrad fields can be viewed as the square root of the metric tensor g^^vix) in the sense of a matrix equation 
[ 2 ^. Inverting the relation Eg. lTMl) we obtain 


V&0 = 90Ax)e%{x)e'<{x). 


(77) 


There are also the following relations 


e&{x)el{x) = di, 

^%{x)eZ{x) = 
e&(x)e^^(x) = g-0, 

e&^(x) = g^^e^(x), 
e^(x)ea!,(x) = gf,^{x). 


(78) 


In terms of the tetrad fields e^{x) and the Dirac matrices 7 “ in the Minkowski spacetime the Dirac matrices 7 ^( 2 ;) 
are defined by 


A covariant derivative we define as 


l^{x) = e^(a;)7“. 


- r^(a;). 


(79) 

(80) 


The spinorial affine connection r^(a;) is defined by 


(^&ix)ep^.^{x), 

where = §(7“7^ — 7 ^ 7 “) and e^^.^(a:) is given in terms of the affine connection r“j^(x) 

4:nix) = e^Jx)-r°‘^^{x)ei{x). 


(81) 


(82) 


In the Einstein gravity the affine connection r“^i,(x) is equal to r“^y(a::) = {°‘^.v} (see Eg. lfTUl) '). Specifying the 
spacetime metric one may transform the Dirac equation Eq. (I72|) into the standard form 




(83) 


where H is the Hamilton operator. For example, for the static metric ds^ = dt^ — W'^ (dr)'^, where V and W are 
spatial functions, one may show that the Hamilton operator is given by (see 


H = 7^^ 


(VE) (VIE)\ 
W ' ^ \ 2E IE / 


(84) 


In the approach to the Einstein-Cartan gravity with torsion, developed above, the Lagrangian of the Dirac field '>p{x) 
with mass m, coupled to the chameleon field through the Jordan metric g^i, = g^^-, is equal to 


C-m — V — 




(85) 
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where the tetrad fields in the Jordan-frame and the Einstein-frame are related by 


= /e“, 

= r 




( 86 ) 


and the covariant derivative is equal to 

= 9^ — ryj(x), 

f = {%.} + S% /-' /,. = {%.} + ^ <)>,.. (87) 

As a result we get 

( 88 ) 

Below we apply the obtained results to the analysis of the electroweak interactions of the neutron and proton. 


B. Electroweak model for neutron and proton, coupled to chameleon field through the Jordan metric 

~ _ r2 

dtJ-iy — J giJ,u 

For the subsequent application of the results, obtained below, to the analysis of the contribution of the chameleon 
field to the radii of the neutron and the proton and to the neutron /3“-decay we defined the electroweak model for 
the following multiplets 



riR, 


(89) 


where '(/'l = ^(1 — 7 ^) V' and 'ipn = ^(1 + 7 ^) ip. Such a model is renormalisable also because of the vanishing of the 
contribution of the Adler-Bell-Jackiw anomalies Qp + Qe = 0, where Qp = +1 and Qe = —1 are the electric charges 
of the proton and electron, measured in the units of the proton charge e [s^ . 

The fermion states Ea. dHiH) have the following electroweak quantum numbers: Nl ■ {Iw = 5 ,= 1), Pfl : {Iw = 
0 , Fu, = 2 ), nji : {1^ = = 0 ) and £l ■ {Iw = 5 ,^^, = — 1 ), : {Iw = 0,Yw = — 2 ), where the third component 

of the weak isospin Iws and weak hypercharge Yw are related by Q = Iw 3 + Fu,/2. In the Einstein-Cartan gravity 
with the torsion field and the chameleon field, coupled to matter field through the Jordan metric g^i, = g^v, the 
Lagrangian of the fermion fields Eg. (15^ . coupled to the vector electroweak boson fields and the Higgs boson field, is 
equal to 


-^ewf 

7^9 




Nl 


fNL f ^i'l'"{x)(^dp-i-g' f 
+ 7pr f~^'il''{x)(dp-ig'f~^ Bp-Fp^ PR + 7riR iji^{x)(dp - Fp^ 


nR 


+ 7 £l f ^il^{x){dp + i^g' f ^Bp-i^gf ^ f ■ Ap - F p 
+ fe7f~^h'"{x)(dp+ig'f~^Bp-Fp)] e^. 


£l 


(90) 


The masses of the neutron and electron one may gain by virtue of the following interactions with the Higgs field $ 




Kn Ur + Ur^'^Nl) - f'^ Re (II^ + C^^'I'^l) 


r4 - r4 _ r4 -- r4 -- 

= -j mnUn-f -^nnip-f mg e e -/ e e (p, 


(91) 
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where k„ and Kg are the input parameters, defining the neutron and electron masses = Knvjy/2 and mg = Kgvj 
respectively. In principle, the proton mass we may gain due to the interaction of the proton fields with the Higgs field 
$, defined by the column matrix with the components {v + ip)(y/2 and Using the Higgs field l> one gets 


6 C 


hp 




-f'^ Kp^Nl<^Pr = -f^nippp- ^ppif 


(92) 


where nip = Kpvjy/2 is the proton mass. In terms of the physical vector field states the Lagrangian of fermion field 
is given by 

1 ^ 

fpix) i-i^[x){dp - i f~^eAp{x) + i 2 cos gw - 2 sin^ 9w)Zp{x) - r^(x)) (— p{x) 


r 

■^cw 

7^9 


+ fp{x)i-i^{x){dp-if ^eAp{x) + if ^ 
+ fn{x)i-i^^{x)(^dp-if 


2 cos 6w 
9 


cos 9w 


2 cos 9w ^ 
•l+yS 


Zpix) - 


sin^ 9w Zp{x) - rp(x)^ p{x) 

1 - 

' n{x) 


+ fn{x)i-f>^{x)(dp-rp,{x)^ (—^)n(x) 

+ f ri{x)W+{x)+f-^n{x)-f^^{x) ^ )p(x)IU~(x) 

1 5 

+ f j>e{x)i-t^{x){dp + i f~^ ^ ^ T- i^\\ f 'y 


2cos9w^^^ ' 


iyg{x) 


Zpix) -Tpix) 

+ fe~{x)i'yf^{x)(^dp+if~^eAp{x) - i f~^ 2 cos 9w ^^ - 2 sin^ 9w) Zp{x) - r^(a;)) ) e~{x) 

+ fe~{x)ij^{x)(^dp+if~^eAp{x) + i f~^ sm'^ 9 w)Zp{x) - r^(a;)) ) e~{x) 


cos 9\y 

- nipp{x)p{x) - m„ n{x)n{x) - nig e~ (x)e“ (x) 

- f^P{x)p{x) p{x) -7^ n{x)n{x) p{x) - 7 ^ e~ {x)e~ {x) p{x). 

We note that for the calculation of r^(x) we have to use the affine connection, given by Eq. dm). 


(93) 


VII. CONTRIBUTION OF CHAMELEON FIELD TO CHARGE RADII OF NEUTRON AND PROTON 

AND TO NEUTRON d"-DECAY 

The electroweak model in the Einstein-Cartan gravity with the torsion and chameleon fields, analysed above, is 
applied to some specific processes of electromagnetic and weak interactions of the neutron and proton to the chameleon 
field. In this section we calculate the amplitudes of the electron-neutron and electron-proton low-energy scattering 
with the chameleon field exchange and define the contributions of the chameleon field to the charge radii of the 
neutron and proton. We calculate also the contribution of the chameleon field to the energy spectra of the neutron 
/3“-decay and the lifetime of the neutron. The calculations we carry out in the Minkowski spacetime replacing metric 
tensor in the Einstein frame by the metric tensor of the Minkowski spacetime ^ giiv The Lagrangian of 

the electromagnetic and electroweak interactions of the neutron, the proton, the electron and the electron neutrino, 
coupled to the torsion field and the chameleon field in the Minkowski spacetime, is given by 

-Cew = fp{x) - i f~^ e Ap{x) +if~^ 4,cL9w ~ 9w)7)Zp{x) - r^(x)) - mp/|p(x) 

+ f n{x) - i f-^ “ 7) Zp{x) - r^(x)) - m„ / | n{x) 

+ f Ux) {d, + ^ J^Zpix) - r,(x)) ( 1 ^) .g{x) 

+ fe~{x) 1*7^ ( 5 ^ +i eAp{x) -if~^ 4 , cos 9w ~ ^sin^gvv) - 7)Zp{x) - r^(x)) - me/| e“(x) 

C 2 ) ^ ^ h(x) 7 ^ ^ ) p{x) W-{x) 
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e' K^+q e' 



FIG. 3: Feynman diagram for the contribution of the torsion-chameleon field to the squared charge radius of the neutron 




1-^5 


V2 


V2 


(94) 


where the ellipsis denotes the interactions with the Higgs field. Then, rup = 938.27205 MeV, = 939.56538 MeV and 
nie = 0.51100MeV are the masses of the proton, neutron and electron, respectively [l^. Expanding the conformal 
factor / = in powers of the chameleon field and keeping only the linear terms we arrive at the following 

interactions 


^int — 


+ 

+ 

+ 


- p{x)p{x) (j){x) — (3n{x)n{x) (j){x) — Pe {x)e (x) cj){x) 


Mpi' 


Mpi 


Mpi 


ep{x)j>^p{x) Ap{x) - 4 cog 5 )^ - (1 - 4: sin^ 9wh^'jp{x)Zp{x) 


ee (a;) 7 ^e (a;)H^(a:) + 
9 

4 cos 9\y 


9 


4 cos 9w 
{x)-f^ (l - j^^n{x)Zp{x) 

'1-7' 


P(a;) 7 '"((l- 4 sin^ 6 »vi/)- 7 '')e {x)Zp{x) 


9 


4 cos 9w 


z/e(a:)7'"(l - l^)ve{x) Zp{x) 
1-75 


-^p(a;) 7 ^ (— n{x)W+{x) + -^ n(a:) 7 ^ (— i^')pi.x)Wp {x) 


V 2 

9 

V 2 ' 


{x)Y ^ ) e (x) W+ {x) + -^e {x)j>^ ) Ve{x) (a;) 


1-75 


V2 


(95) 


where we have omitted the interactions with the Higgs field, which do not contribute to the processes our interest. 
The contribution of the terms, containing r^(x), which in the Minkowski spacetime is equal to T^ = | (T^i^(£n/)’‘^, 
can be transformed into total divergences and omitted. 


A. Contributions of chameleon to squared charge radius of neutron 


The torsion (chameleon) contribution to the squared charge radius of the neutron is defined by the Feynman diagram 
in Fig. [21 The chameleon-neutron interaction is defined by (see Ea. dMll 'l 

Tin 

Cnn<i>{x) =-j3-^n{x)n(x)(j){x). (96) 

Mpi 

The Lagrangian of the interaction, given by Eq. o for \/—g = 1 , can be transcribed into the form 


C^j^{x) = -I- {x)Fp^{x)<j){x), 


(97) 


2 Afpi 

where we have omitted the total divergence. The analytical expression for the Feynman diagram in Fig. [3] is equal to 

P^TTln f d'^q ,, 1 


M(e“ 


n —>■ e n 


) = 


M|i J {2nyi L 




TOe — ke — q — iO 


Yu{ke) (?7"\^-7"‘"9^)44..(g) 


^ {Vax{q T ke ke}tp 9atp(^q F ke t^{q F ke fcg) 


1 


ml - [kn - k'^y 


[u{k'Ju{kn)], 


(98) 


where is the chameleon mass Ea. (l39|) as a function of the chameleon-matter coupling constant /?, the environment 
density p and the Ratra-Peebles index n, and Dap{Q) is the photon propagator Ea. (|i5)) . Substituting the photon 
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propagators D°^^{q) and Dav{q + ke — /c'), taken in the form of Ea. ipS)) . into Ea. dMl) one may show that the integrand 
does not depend on the gauge parameter i.e. the integrand is gauge invariant. 

Measuring the electric charge on the neutron in the electric charge of the proton for the calculation of the neutron 
electric radius we have to compare Eg- dMl) to the amplitude 


M(e n->• e n) =[w(fcg)u(fce)][M(A:^)w(fc„)], (99) 

where e, = —Cp and Cp are the electric charges of the electron and proton, respectively, and is the squared charge 
radius of the neutron. The product [u{k'g)u{ke)][u{k'^)u{kn)] is equivalent to the product [u{k'^)^^u{ke)]{u(k'^)^'^u(kn)] 
in the low-energy limit. 

From the comparison of Ea. dMl) with Ea. dMl) the contribution of the chameleon to the squared charge radius of 
the neutron can be determined by the following analytical expression 


2 6 / 3^771 r d^q 

J {2Tr)H 

where we have set fc' = ke and k'^ = k^. 


l{ke)^^ 


me — ke — q — iO 


ru{ke) iv^'^q^-q^ql 


{q^ + * 0 ) 


Merging denominators by using the Feynman formula 

1 2xdx 

^ Jo [Ax + B{1 - x)]3 


2 ’ 


we arrive at the following expression 

6/3^to 


= 




dx 2x 


d^q u{ke)j^ime + ke + q)Yu{ke) 


(2TT)'^i [mg(l — x)'^ — id + keil — 




( 100 ) 


( 101 ) 


( 102 ) 


Making use a standard procedure for the calculation of the integrals Ea. H102p . i.e. i) the shift of the virtual momentum 
q+keil — x), ii) the integration over the 4-dimensional solid angle and iii) the Wick rotation, we arrive at the expression 


2 6/3^mem , r, f 9(1 — a;)q^ — 3mg(l — cc)^ 9 P^mem ( ^ \ 

■ mlM^, Jo J W [ 92 +m 2 (l-x) 2]3 = " 4^^ t j ’ 

where M is the ultra-violet cut-off. For numerical estimates we set M = Mpi [2^. This gives 

2 _ 9 /3^ mem / Mpi \ 

” 477^ Mpj TOg / 

According to [s^ , the squared charge radius of the neutron can be defined by the expression 


ma ’ 


(103) 


(104) 


(105) 


where a = 1/137.036 and 6„g are the fine-structure constant and the electron-neutron scattering length, respectively. 
For the experimental values of the electron-neutron scattering lengths 6„g = (—1.330 ± 0.027 ± 0.030) x 10“^ fm and 
hne = (—1-440 ± 0.033 ± 0.030) x 10“^ fm, measured from the scattering of low-ei ierg y electrons by ^°®Pb and ^°®Bi 
[ 33 , respectively, we get (r^)exp = —0.115(4) fm^ and (r^)exp = —0.124(4) fm^ [sj, respectively. The theoretical 
value of the squared charge radius of the neutron is 


«2 

rl = -3.588 X 10"^^ ^ W = -6.229 x 10"^® 
mi 


/3n+l 


/ nAIp\h? \ n+i 


n(n + 1) V 


fm^ 


(106) 


where the chameleon mass is measured in meV. From the comparison to the experimental values we obtain 


/3 > [1.864 X 10^® 
P > [1.991 X 10^® 


n+1 

n(n +1)] n 1 

( Pm, ^ 

n+2 

1 " 

for 

^ nMp\A^ ) 

I 

n+2 

n+1 

n(n + 1)] n 1 


( Pm ^ 

1 " 

for 

^ nMp\A^ J 

1 


208 


Pb, 


208 


Bi, 


(107) 
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FIG. 4: The lower bound of the chameleon-matter coupling constant /? from the experiments on the electron-neutron scattering 
length for the scattering of the slow neutron by lead (left) and bismuth (right), respectively. The shaded area is excluded. 


respectively. One may use Ea. (ll07|) for the estimate of the lower bound of the chameleon-matter coupling constant 
p. Since the experiments on the measuring of the electron-neutron scattering length have been carried out for liquid 
lead and bismuth [s^ with densities ppb = 10.678 g/cm^ and psi = 10.022 g/cm^, respectively, in Fig.|4]we plot the 
lower bound of the chameleon-matter coupling /3 at which the contributions of the chameleon field are essential. 

The minimal lower bound /3 > 10^® is ten orders of magnitude larger compared to the value P < 5.8 x 10®, measured 
recently in the qBounce experiments [l^ . As a result, the contribution of the chameleon field to the electron-scattering 
length in the environment of the liquid lead and bismuth is negligible. Thus, in order to obtain a tangible contribution 
of the chameleon to the electron-neutron scattering length bm or the squared charge radius of the neutron r^, the 
experiments should be carried out in the environments with densities of order p ^ 10“^g/cm® or even smaller. 


B. Contributions of the chameleon field to the squared charge radius of the proton 


The results obtained above for the squared charge radius of the neutron can be applied to the analysis of the 
contribution of the chameleon (torsion) to the squared charge radius of the proton. Since the interaction of the 
chameleon field with the proton is described by the Lagrangian Eg. (1961) with the replacement n{x) —>■ p(a;), where 
p{x) is the operator of the proton field, the contribution of the chameleon field to the squared charge radius of the 
proton Tp is dehned by Eas. ((10411 and (14511 with the replacement —?► rUn rup and rUe —t mp, where nip is the 

mass of the -meson [T3 |. 

The contribution of the chameleon field to the charge radius of the proton has been recently investigated by Brax 
and Burrage . According to Brax and Burrage [s^ , the contribution of the chameleon field may solve the so-called 
“the proton radius anomaly” [sol-li^. As has been shown in the Lamb shift AE 2 s^ 2 p of the muonic hydrogen, 

calculated in QED with the account for the nuclear effects, can be expressed in terms of the charge radius of the 
proton Tp 

AE2s^2p = 209.9779(49) - 5.2262 + 0.0347 r®, (108) 

where AE 2 s-^ 2 p and Vp are measured in meV and fm, respectively. The charge radius of the proton rp = 0.8768(69) fm, 
measured from the electronic hydrogen [d^ and agre eing well with the charge radius of the proton Vp = 0.879(8) fm, 
extracted from in the electron scattering experiments |49l| . In turn, the experimental value of the charge proton radius 
Tp = 0.84087(39) fm, extracted from the measurements of the Lamb shift of muonic hydroge n [ ^ . is of about 4% 
smaller compared to the charge radius of the proton, measured from the electronic hydrogen |4nl |. The correction to 
the Lamb shift, caused by the correction to the squared charge radius Sr^, is equal to 

SE2s^2p = (-5.2262 + 0.0521 rp) Sr^ = -5.180 (109) 

where we have set Vp = 0.8768fm. The correction to the charge radius of the proton, caused by low energy pp 
scattering, is equal to (see Ea. (ll04|) ) 


6 r 


2 

p 


9 P'^ rUpiTip / Mpi \ 

47r2 M|j ^\mp) 


= -6.617 X 10"^® 




( 110 ) 


where = 105.6584MeV and nip = 938.2720MeV are the muon and proton masses, respectively [T^, and mp and 
Svp are measured in meV and fm^, respectively. Substituting Ea. dllOl) into Eg. (11091) we express the correction to the 
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FIG. 5: The chameleon-matter coupling constant /3 as a function of the index n, fitting the experimental value 5 E2s-2p = 
Ai 52 s- 2 p —Ai 52 s- 2 p|r-p=o .8768 = 0.311 meV, where AE 2 s- 2 p = 206.2949 meV [s^, of the Lamb shift of muonic hydrogen [3l.The 
shaded area is excluded. 


Lamb shift of the muonic hydrogen in terms of the parameters of the chameleon field theory and the matter density 
Pm- We get 
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( 111 ) 


with dE 2 s^ 2 p = 0.311 meV [i^. In Fig. [^we plot the coupling constant ^ as a function of the Ratra-Peebles index n 
at the environment density pm ~ lOg/cm^ p. One may see that for SE 2 s^ 2 p = 0.311 meV the lower bound of the 
chameleon-matter coupling constant /3, at which the contribution of the chameleon is tangible, is B > 10^®. This is 
seven orders of magnitude larger compared to the recent experimental upper bound j3 < 5.8 x 10® [12j| . 


C. Contribution of the chameleon field to the neutron /3 —decay 

In this section we investigate the neutron /3“-decay, caused by the interaction with the chameleon. This means 
that we investigate two reactions: i) the neutron /3“-decay with an emission of the chameleon n—>'P + e“+i>e + 0 
and ii) the chameleon induced neutron /3“-decay </> + n —>■ p + e“ + i/e- Since formally these two reactions are related 
by ^ where is a 4-momentum of the chameleon, we give the calculation of the amplitude of the neutron 
/3“-decay n—>'P + e“ + i/e + </' with an emission of the chameleon. 


Neutron fi —decay with the chameleon particle in the final state n ^ p + e -^Vs + cj) 

The calculation of the amplitude of such a decay we use the following effective interactions 

/leff = -p ^^^n{x)n{x)(t){x) - P ^^^p{x)p{x)(j){x) 

Kd |[p(x)7^(l + \-t^)n{x)] + -^d''[p{x)a^,pn{x)]^[e-{x)-i^{l - (112) 

where Gp = 1.1664 x 10“^ MeV“^ is the Fermi coupling constant, Vud = 0.97427(15) is the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix element [l^, A = —1.2750(9) is the axial coupling constant [s^ (see alsof^l 
and K, = K,p — Kn = 3.7058 is the isovector anomalous magnetic moment of the nucleon, defined by the anomalous 
magnetic moments of the proton Kp = 1.7928 and the neutron = —1.9130 and measured in nuclear magneton [l3l| . 
M = (m„ + mp)/2 is the nucleon average mass. 

The Feynman diagrams of the amplitude of the neutron /3“-decay n^p + e~+Df. + ([> are shown in Fig.[6l For 
the calculation of the analytical expression of the decay amplitude we follow (s^ and carry out it in the rest frame 
of the neutron, keeping the contributions of the terms to order 1/M. The energy spectrum and angular distribution 
of the neutron /3“-decay with polarised neutron and unpolarised proton and electron we may write in the following 
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FIG. 6: Feynman diagrams of the neutron fi -decay with an emission of the chameleon n^p + e -\- Ue + cf). 


general form 


dri 2 


Gl 


pe- i^e F{Ee,Z = 1) ^(ke,K,k^) - kp - ke - k^ - k^), (113) 


4m„ 


Mpj 


where F{Ee,Z = 1) is the relativistic Fermi function, describing the final-state electron-proton Coulomb interaction 
[a, is the phase-volume of the decay final state 


d^kp d^ke d?ku d?k^ 
{2nf2Ep {2-Kf2Ee {2-Kf2E^ {2 ttY2E^ 


(114) 


and kj for j = and ^ is a 4-momentum of the neutron and the decay particles, respectively. The factor 

^(ke,kv,kff,) is the contribution of the phase-volume, taking into account the terms of order 1/M. Following one 
obtains 


4>(fce, k^, k^) = 1 + ^ ^Ee 



k(p ' ki/ \ 

E, ) 


2 E^E(j) ke * kfj) 
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(115) 


The last two terms define the deviation from the phase-volume factor, calculated in (s^ for the neutron (3 -decay 
n —^ p + e“ + f'e- Taking into account the phase-volume factor 4>(fce, k^) we may carry out the integration over 

the phase-volume of the n—>'P + e“+i/e + ^ decay, neglecting the contribution of the kinetic energy of the proton. 

Then, |M(n pe~ />)P is the squared absolute value of the decay amplitude, summed over the polarisation of 
the decay electron and proton. The analytical expression of the amplitude M{n ^ pe~ </) is defined by 


M{n ^ pe r'e />) 
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'Up(/up,(Jp) A ,v ^ n i^k f, ^ 
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rUn - kn + k^ - iO J 

L z J 


(116) 


where Uj{kj, Uj) for j = n,p and e are the Dirac bispinors of fermions with polarisations ctj, {k^, +1/2) is the Dirac 
bispinor of the electron antineutrino and Op is defined by [1^ 


Ol^\kp, kn) = 7^(1 + Ay*^) + i — ap^{kp - ± k^Y. 

In the accepted approximation the amplitude Eq. (| 115|) can be defined by the expression 

2 r 1 
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where is given by 


Of, = 7p(l + Xy^) + i — af,^{kp - k^f 
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and the matrices Oq and O, taken to order 1/M, are equal to 

A 

i-^ + - 

= 


-1-W + 2M<"’‘A) 


and 


\ 2M J 2M 2M ^ 

a 



\ 


—Act 

:-S)- 

i^+i7^{dxk, 
2M 2M ^ ^ 

b 


( 120 ) 


O = 


( 121 ) 

where Eq = ((m„ — m^)^ — + mg)/2(m„ — rricj)) is the end-point energy of the electron-energy spectrum. The 

matrices 0° and O are defined to order 1/M only . For the calculation of the amplitude of the /3“-decay of the 
neutron we use the Dirac bispinorial wave functions of the neutron and the proton 


i(0^ (Jn) — \/2y7ryj ^ Q ^ ’ '^pi ; ^p ) — \/F'p T ' 


^P 

a ■ kp 

Ep + ruj 


'fip 


( 122 ) 


where ifn and (pp are the Pauli spinor functions of the neutron and proton, respectively. For the energy spectrum and 
angular distribution of the neutron /!“-decay with polarised neutron and unpolarised proton and electron we may 
write in the following general form 
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(123) 


where ke = -^ZFi/ — rUg is the absolute value of the electron 3-momentum and Zn is the unit polarisation vector of 
the neutron |Zra| = 1- The correlation coefficients ({Eg), a{Ee), A{Ee), B{Ee), Kn{Ee), Qn{Ee) and D{Ee) can taken 
from [ 5 ^ at the neglect of the radiative corrections. The correction coefficient we may represent in the following 
form 


Fa — F^^^ 


. p(3) 

^4> ' -^ 4 , ^ 


(124) 


where the correlation coefficients E^^'^ for j = 1,2,3 are defined by i) the contributions of the terms, depending on 
the energy and 3-momentum of the chameleon particle in the matrices Oq and O given by Eos. (11191) and (jl20L 
ii) the dependence of the 3-momentum of the proton kp on the 3-momentum of the chameleon particle, caused by 
the 3-momentum conservation kp = —ke — k^ — k^, (see Eqs.(A.16) and (A.17) in Appendix A of Ref.jl^) and hi) 
the contributions of the phase-volume factor Eq. (HHl), respectively. The analytical expressions of these correlation 
coefficients are equal to 
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and 
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(127) 


Now we may integrate over the phase-volume of the n ^ p + e + Vp + (f) decay. First of all we integrate over the 
3-momentum of the proton kp. As has been mention above, we make such an integration at the neglect of the kinetic 
energy of the proton. Then, since one can hardly observe the dependence of the energy spectrum and the angular 
distribution on the direction of the 3-momentum of the chameleon, we make the integration over the 3-momentum 
of the chameleon k^. The obtained energy spectrum and angular distribution is 

d^Xn^iEp, E^, E^, kp, fc., Cn) _ ^ ^ -Ep-E^- E^) kpEp El (1 + 3A2) C^{Ep) 
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The correlation coefficients ((j,{Ep), a^{Ep), A^{Ep) and B^{Ep) are equal to 

Epk 


C4Ee) = aEe) + 


M'' 


a^{Ep) — a{Ep) + ao (^1 + 2 ^ 

A(p{Ep) = A{Ep) + (^-2X +Aq (^1 + 2 
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(129) 


where oq = (1 — A^)/(l + 3A^), Aq = —2A(1 + A)/(l + 3A^) and Bq = —2A(1 — A)/(l + 3A^) [5^(see also Hi). The 
correlation coefficients C{Ep), d{Ep), A{Ep) and B{Ep) are calculated in by taking into account the contributions 
of the weak magnetism and the proton recoil to order 1/M but without radiative corrections. 

The rate of the decay n^p + e~ + Dp+ 4> diverges logarithmically at ^ 0. We regularise the logarithmically 
divergent integral by the chameleon mass m^. As result we get 


27r3 


(130) 


where /^{Eq, Z = 1) is the Fermi integral 


U{Eo,Z =1) = f ^dEp Ep ^El - ml (Aq - Epf \in( ^° ^ i ^Ee) E{Ep, Z = 1), (131) 
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where C{Ep) is given by (see Eq.© of Ref.js^) 
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In Fig. [7] we plot the electron-energy spectrum p^{Ep) of the neutron /3 -decay with an emission of the chameleon, 
defined by 
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where f{Eo, Z = 1) is the Fermi integral calculated in [^, and compare it with the electron-energy spectrum p^- {Ep) 
of the neutron /3“-decay calculated in (see Eq.(D-59) of Ref.jl^). The chameleon mass is determined at the 
local density p ~ 1.19 x lO^^i g/cm^ = 5.12 x lO”^'^ MeV"^. This is the density of air at room temperature and pressure 
P ~ 10“3 mbar [^|^. For /3 < 5.8 x 10® we obtain Xn<p < 2 x 10“®^ s”^ and the branching ratio BR „0 < 1.8 x lO”®^ 
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FIG. 7: (left) The values of the Fermi integral f^{Eo, Z = 1) as a function of the index n for n = 1, 2,..., 10. (right) The 
electron-energy spectra p^{E^) (continuous lines) and p^-{E^) (dashed line) of the neutron /?“-decay with and without an 
emission of a chameleon, respectively. The densities p^{Ei,) depend slightly on the index n and are represented by only one 
continuous blue line. 




(a) (b) 

FIG. 8: Feynman diagrams for the reaction (j> + n ^ p + e~ + Ue 


Neutron {3 —decay (j) -\- n ^ p -\- e v^, induced by the chameleon field 


For the calculation of the amplitude M{(j)n ^ pe~ 9^) of the induced neutron /3“-decay 4> + n^p + e~ + 9e we 
may use the amplitude Ea. (lll6|) with the replacement —>• —where is a 4-momentum of the chameleon. The 
Feynman diagrams for the chameleon-induced neutron /3“-decay is shown in Fig.|Sl The cross section for the induced 
neutron /3“-decay is 
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where the contribution of the electron-proton final-state Coulomb interaction is not important and neglected. Using 
the results, obtained in previous subsection, we get 
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Integrating over we obtain 
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Using the results, obtained in (see also [b^), we may analyse the quantity 
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which defines the number of transitions n ^ pe Vg per second, induced by the chameleon, where $ch(£' 0 ) is the 
number of solar chameleons per eV • s • c m^, normalised to 10% of the solar luminosity per unite area Lq/AttRq = 
3.9305 X lO^^eVs-i cm’^ = O.OlOOTeV^ [gi, where Lq = 2.3893 x lO^^gys"! and Rq = 6.9551(4) x 10^°cm are the 
total luminosity and the radius of the Sun |l3j|. Following we obtain that X^n < 10“"^^ s“^ for /3 < 5.8 x 10® [T^ . 


VIII. CONCLUSION 

We have developed the results, obtained in [l^, where there was shown that the chameleon field can serve also as 
a source for a torsion field and low-energy torsion-neutron interactions, where the torsion field is determined by a 
gradient of the chameleon one. Following Hojman et al. [l^ we have extended the Einstein gravitational theory with 
the chameleon field to a version of the Einstein-Cartan gravitational one with a torsion field. For the inclusion of the 
torsion field we have used a modified form of local gauge invariance in the Weinberg-Salam electroweak model with 
minimal coupling and derived the Lagrangians of the electroweak and gravitational interactions with the chameleon 
(torsion) field. 

Gauge invariance of the torsion-photon interactions has been explicitly checked by calculating the amplitudes of 
the two-photon decay of the torsion (chameleon) field 7 + 7 and the photon-torsion (chameleon) scattering 

7 + (()—^0 + 7 or the Compton photon-torsion (chameleon) scattering. Unlike the Compton-scattering, where 
photons scatter by free charged particles with charged particles in the virtual intermediate states, in the photon- 
torsion (chameleon) scattering a transition from an initial (7 (p) state to a final (7 (p) goes through the one-virtual 
photon exchange (see Fig. [5^ and Fig.j^b) and the local interaction (see Fig.[5J:). The Feynman diagram Fig.j^Ji 

is self-gauge invariant due to the local £ 7^0 interaction. Gauge invariance has been checked directly by a replacement 
of the one of the polarisation vectors of the photons in the initial and final state by its 4-momentum. Since in these 
reactions the coupling constant of the photon-torsion (chameleon) interaction is ^eff = P/Mp\, in analogy with gauge 
invariance of photon-charge particles interactions, where electric charge is a coupling constant - unrenormalisable 
by any interactions, one may assert that the coupling constant ges = P/Mpi should be also unrenormalisable by 
any interactions. This ma y p lace some strict constraints on possible mechanisms of the chameleon-matter coupling 
constant /3/Mpi screening [^[^. In this connection the Vainstein mechanism, leading the screening of the coupling 
constant j3/Mp\ by the factor where Z > 1 is a finite renormalisation constant of the wave function of the 

chameleon field caused by a self-interaction of the chameleon field or some new higher derivative terms [^ . is 
prohibited in such a version of the Einstein-Cartan gravity with the chameleon field and torsion. Because of the 
smallness of the constant /3^/Mpj < 10“®°barn/eV^, estimated for P < 5.8 x 10® [l^, the cross section for the 
photon-chameleon scattering is extremely small and hardly may play any important cosmological role at low energies, 
for example, for a formation of the cosmological microwave background and so on. However, since in our approach 
the coupling constant I3^/Mp\ is fixed in terms of the coupling constant 13/Mp\, the recent measurement of the upper 
bound [3 < 5.8 x 10® can make new constraints on the photon-chameleon oscillations in the magnetic field of the 
laboratory search for the chameleon field [gtI - I^ . 

In our approach the effective chameleon-photon coupling ^eff = P/Mpi is equal to geff = P/Mpi < 2.4x 10“^° GeV“^, 
where we have used the experimental upper bound of the chameleon-matter coupling constant j3 < 5.8 x 10® [l^ . 
The obtained upper bound ^eff < 2.4 x 10“^° GeV“^ is in qualitative agreement with the upper bounds, estimated 
by Davis, Schelpe and Shaw [ 5 ^. Then, the constraints on /3: j3 < 1.9 x 10^ (n = 1), /3 < 5.8 x 10^ {n = 2), 
/3 < 2.0 X 10® (n = 3) and j3 < 4.8 x 10® (n = 4), measured recently by H. Lemmel et al. [s^ using the neutron 
interferometer, may be used for more strict constraints on the astrophysical sources of chameleons, investigated in 

[Hl-p. 

Using the photon-torsion (chameleon) interaction we have estimated the contributions of the chameleon field to the 
charged radii of the neutron and proton. All tangible contributions can appear only for /3 ^ 10®. This, of course, is not 
compatible with recent experimental data j3 < 5.8 x 10® by Jenke et al. [l^. The branching ratio for the production of 
the chameleon in the neutron /3“-decay n —>■ p+e~ +i^e+4> is extremely small Br(n p e~ Dg cp) < 1.8 x 10“®^. In turn, 
the half-life of the neutron = in2l\^n, caused by the chameleon induced neutron /3“-decay p + n —?► p+e~+Vg, 

is extremely large = In^jX^n > 2 x 10®® yr. Of course, because of the neutron life-time = 880.3(1.1) s 

being in agreement with the recent theoretical value = 879.6(1.1) s [5^, the chameleon induced neutron /3“-decay 
cannot be observed by a free neutron. The experiment, which can give any meaningful result, can be organised in a 
way, which is used for the detection of the neutrinoless double /3 decays . For example, it is known that the isotope 
^®Ge is both stable with respect to non-exotic weak, electromagnetic and nuclear decays and are neutron-rich. It 
is unstable only with respect to the neutrinoless double /3“-decay 0 - The experimental anal ysis of the low-bound 
on the half-life of ^®Ge has been carried out by the GERDA Collaboration Agostini et al. [0)[23 measuring the 
energy spectrum of the electrons. The experimental lower bound has been found to be equal to T 1/2 > 3 x 10^® yr 
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(90 % C.L.). We would like also to mention the experiments on the proton decays, carried out by Super-Kamiokande 
For specific modes of the proton decay, i.e. p —>■ e+Tr^, p —>■ /x+7r° and p —>■ for the half-life of the proton 

there have been found the following lower bounds: T 1/2 > 5.7 x 10^^ yr, T 1/2 > 4.6 x 10^^ yr and T 1/2 > 4.1 x 10^^ yr, 
respectively. 

Since the lifetime of the neutron r„ = 880.3(1.1) s [l^, one cannot analyse experimentally the chameleon-induced 
/3“-decay on free neutrons. For the experimental investigation of the chameleon-induced /3“-decay one may propose 
the following reaction 

(j> + ^ 4 gCd —>■ ^ 4 gln + e + be, (138) 

where ^IgCd is an atom with a stable nucleus in the ground state with spin^P™*^^ J'^ = O’*'. Since ^Igln is an atom 
with a nucleus in the ground state with spin^P®'*''*^) = 1+, the reaction Eg. (113811 is determined by the Gamow-Teller 

transition 'i|Cd —>■ ^l^ln for chameleon energies > 3.088 MeV. The calculation of the threshold energies in weak 
decay of heavy atoms with the account for the contribution of the electron shells can be found in . 

We have to note that the atom ^Ifln is unstable under the electron capture (EC) and /?“ decays with the branches 
56% and 44%, respectively [t^. Since the EC decay of ^ifln is not observable in the experiment on the chameleon- 
induced /3“-decay, one may observe the /?“~decay of ^^In. However, the electron energy spectrum of ^ifln is restricted 
by the end-point energy Eq = 0.673 MeV and can be distinguished from the energy spectrum of the electron, appearing 
in the final state of the reaction Eg. (11381) . 

The main background for the chameleon-induced /3“-decay Ea. (ll38p is the reaction 

I'e + ^4gCd —> ^Igln + e , (139) 

caused by solar neutrinos with energies > 3.088 MeV. Because of the threshold energy the reaction Eg. (11391) can 
be induced by only the solar ®B and hep neutrinos (l3|. S ince the /lep-solar neutrino flux is approximately 1000 times 
weaker in comparison to the ®B-solar neutrino flux |l3l| . the reaction Eo. (11391) should be induced by the ®B-solar 
neutrinos. 

Concluding our analysis of standard electroweak interactions in the gravitational theory we would like to discuss 
the results, obtained recently by Obukhov et al. [t^. There, the behaviour of the Dirac fermions in the Poincare 
gauge gravitational field including a torsion was analysed. The Hamilton operator of the spin-torsion interaction has 
been derived [^ . In a weak gravitational field and torsion field approximation, which we develop in this paper, such 
a spin-torsion interaction takes the form 

Hspi„-tors = -i(S-f Ty^T'’), (140) 

where E = and y® = zy^y^y^y^ are the Dirac matrices [ 2 ^. Then, and T are the time and spatial 

components of the axial torsion vector field T°‘ = (T°, T ), defined by 

= (141) 

where T/ 3 fj,i, is the torsion tensor field and is the totally antisymmetric Levi-Civita tensor = 1 [^. Using 

the experimental data and on the measurements of the ratio of the nuclear spin-precession frequencies of 
the pairs of atoms (^®®Hg, ^°^Hg) and (^He, ^^®Xe) [t^ with nuclear spins and parities (= 4 , J'^ = | ) and 
(J^ = 4 , = 4 )^ respectively, Obukhov, Silenko and Teryaev have found the strong new upper bound on the 

absolute value of the torsion axial vector field T. They have got 

ifllcosei < 4.7 X 10"2^eV. (142) 

In the approach, developed in our paper, the tensor torsion field is equal to = {13/M-p\){gp^d^4> — gp^dv(j>) 

(see Eq.([71)). Multiplying such a tensor torsion field by the totally antisymmetric Levi-Civita tensor gg|; 

T“ = 0. Thus, the upper bound of the absolute value of the axial vector torsion field Eq. (HH), obtained by Obukhov, 
Silenko and Teryaev [73|, does not rule out a possibility for a torsion field to be induced by the chameleon field as it 
is proposed in our paper. 
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